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Abstract 

The decay J/V' ^ NNtt provides an effective isospin 1/2 filter for tfie ttN system due to 
isospin conservation. Using 58 million J/ijj decays collected with the Beijing Electromagnetic 
Spectrometer (BES) at the Beijing Electron Positron Collider (BEPC), more than 100 thousand 
J/if) ^ pTr~n + c.c. events are obtained. Besides the two well known N* peaks at around 1500 
MeV/c^ and 1670 MeV/c^, there are two new, clear N* peaks in the pTr invariant mass spectrum 
around 1360 MeV/c^ and 2030 MeV/c^ with statistical significance of llcr and 13cr, respectively. 
We identify these as the first direct observation of the A/'*(1440) peak and a long-sought "missing" 
A''* peak above 2 GeV/c^ in the ttN invariant mass spectrum. 

PACS numbers: 14.20.Gk, 13.75.Gx, 13.25.Gv 
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The nucleon is the simplest system in which the three colors of QCD can combine to form 
a colorless object, and the essential nonabelian character of QCD is manifest. It is necessary 
to understand the internal quark-gluon structure of the nucleon and its excited N* states 
before we can claim to really understand the strong interaction. 

A very important source of information for the nucleon internal structure is the A^* mass 
spectrum as well as its various production and decay rates. Our present knowledge of these 
comes almost entirely from irN experiments performed more than twenty years ago Q|. 
Because of its importance for the understanding of nonperturbative QCD, a series of new 
experiments on A^* physics with electromagnetic probes (real photons and electrons with 
space-like virtual photons) have recently been started at JLAB, ELSA at Bonn, GRAAL at 
Grenoble and SPRINGS at JASRI 2]. They have already produced some important results 
Nevertheless, our knowledge of the A^* resonances remains very poor. Even 
for the well-established lowest excited state, the A^*(1440), properties such as mass, width 
and decay branching ratios etc., still have large experimental uncertainties Q]. Another 
outstanding problem is that, in many of its forms, the quark model predicts a substantial 
number of "missing" A^* states around 2 GeV/c^, which have not so far been observed 
The difficulty in extracting information on these high mass A^* resonances from vrA^ and 7A^ 
experiments is the overlap of many broad resonances with various spin and isospin. Recently 
J/ip decays at the BEPC were proposed as an excellent place to study A^* resonances, 
and A^* production from J/tp ^ ppr] has been studied based on 7.8 million BESI J/tp events 

In this letter, we report on a study of A^* resonances from J/ip —>■ pir^n and pTT~^n 
channels based on 58 million J/ip events collected with the BESII detector at the BEPC. 
Due to isospin conservation, the vrA^ system in J/tp decay is pure isospin 1/2. Compared 
with vrA^ and 7A^ experiments, which mix isospin 1/2 and 3/2, this is a big advantage. 

BESII is a large solid-angle magnetic spectrometer that is described in detail in Ref . [l^ . 
Charged particle momenta are determined with a resolution of ap/p = 1.78%a/1 + p"^ (withp 
in GeV/c) in a 40-layer cylindrical drift chamber. Particle identification is accomplished us- 
ing specific ionization {dE/dx) measurements in the drift chamber and time-of-fiight (TOP) 
measurements in a barrel-like array of 48 scintillation counters. The dE/dx resolution is 
o'dE/dx = 8.0%, and the TOP resolution is o"tof = 180 ps for the Bhabha events. The 
combination of dE/dx and TOP information provides efficient identification of the charged 
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pion and proton (or antiproton) for the J /if) ^ p7T~n and p-K^n processes. 

A GEANT3 based Monte Carlo simulation package (SIMBES)^^ with detailed consid- 
eration of real detector performance (such as dead electronic channels) is used for estimating 
the detection efficiency and background contributions. The consistency between data and 
Monte Carlo has been carefully checked in many high purity physics channels, and the 
agreement is quite reasonable. 

For the decay J/ip pTT~n, the anti-neutron is not detected directly. We select the p 
and 7r~ from two prong events with oppositely charged tracks and require the missing mass 
to be consistent with the n mass. 

In the event selection, we require each charged track to be well fitted to a helix originating 
near the interaction point and to be within the polar angle region | cos 6*1 < 0.8. Both the 
dE/dx and TOF information are used to form particle identification confidence levels p^^^, 
where i denotes vr, or p. For the positive charged track, we require p^^^ > p^^^ and 
Plid > Ppid, for the negative charged track, we require p^^^ > Ppid and p^^^^ > p^^^. After 
these requirements, the p7T~ invariant mass spectrum shown in Fig.^is obtained. There is a 
narrow A peak at 1.116 GeV due to background channels with As. After imposing a further 
requirement of Mp^ > 1.15GeV/c^ to remove this background, the missing mass distribution 
shown in Fig. |21 is obtained. A clear n peak is observed. 




Mp^-(GeV/c^) 

FIG. 1: The pn invariant mass spectrum before the Mp,r > 1.15GeV/c^ requirement. 

A similar selection is used for the charge conjugate channel J/ip ^ pir^n. The missing 
mass distribution of pvT"*" is almost identical with that for J /iIj ^ pir'n. Fitting the missing 
mass spectra between 0.60 and 1.15 GeV/c^ with signal plus background described by a 
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FIG. 2: The missing mass distribution for J/il^ pir n with the fitted background shown by the 
shaded area. 

fourth order Chebyshev polynomial and using Monte-Carlo determined acceptances for these 
two channels, branching ratios for these two channels are measured as (2.36 ± 0.02 ± 0.21) x 
10"^ and (2.47 ± 0.02 ± 0.24) x 10~^ for p7r~n and pvr+n, respectively. The statistical 
errors are small due to the large data sample. Systematic errors come from several sources. 
Background uncertainty is about 4%. The difference of particle identification between data 
and Monte Carlo is about 5%. Uncertainty on the total number of J/ip events is ~ 4.7%, and 
the tracking error is ~ 4%. Combining these errors in quadrature gives a total systematic 
error of (9 - 10)%. 

Taking n/n events in the region of the missing mass distribution within 0.94 ± 0.06 
GeV/c^, we get 58822 and 54524 events for J/ip p7c~n and p7c~^n, respectively, including 
about 6% background. The Dalitz plots for these two channels are shown in Fig. El and are 
very similar. The asymmetry between pvr and mr is partly due to differences in detection 
efficiency and partly due to isospin symmetry breaking effects from the electromagnetic 
interaction. The pvr and nir selection efficiencies versus invariant mass in the decay J/ip ^ 
pTT~n are shown in Fig. |3 The efficiency has a sharp drop for ravr" invariant mass above 
2 GeV/c^, where the recoil proton has low momentum and BESII detection efficiency is 
correspondingly low. 

Shown in Fig. El are the pvr" and pvT"*" invariant mass spectra, as well as the phase space 
distributions (times detection efficiency) obtained with the SIMBES Monte Carlo. The two 
invariant mass spectra look very similar. 
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FIG. 3: Dalitz plots of M^^ vs. for J/^; pTT~n (left) and pvr+n (right). 
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FIG. 4: The selection efficiency versus invariant mass for pir (upper plot) and mr (lower plot) in 
the decay J /il) ^ pTr~n . 

To investigate the amplitude squared behavior as a function of invariant mass, we remove 
the phase space factor and efficiency factor from the invariant mass distribution by dividing 
the data by Monte Carlo phase space times the detection efficiency. The results are shown 
in Fig. ini with data and Monte Carlo for the two channels normalized to the same number of 
events as for pir^n data. The results for the two charge conjugate channels are consistent. 
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FIG. 5: The pir and pn^ invariant mass spectra for J/i/; — > pvr n (left) and pn^n (right), 
compared with phase space. 

At low pir invariant mass, the tail from the nucleon pole term, expected from theoretical 



considerations 
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, 1^, is clearly seen. There are clearly four peaks around 1360, 1500, 1670 
and 2065 MeV/c^. Note that the well known first resonance peak (A(1232)) in vrA^ and 
scattering data does not show up here due to the J/ip decay isospin filter. While the two 
peaks around 1500 MeV/c^ and 1670 MeV/c^ correspond to the well known second resonance 
peak and the third resonance peak observed in vrA^ and scattering data, the two peaks 
around 1360 MeV/c^ and 2065 MeV/c^ have never been observed before in vrA^ invariant 
mass spectra, although a peak around 1400 MeV was observed in the missing mass spectrum 
in pp ^ p + X reaction |l4| with X containing not only vrA^ but also ttttN, StcN, etc.. The 
one around 1360 MeV/c^ should be from A^*(1440) which has a pole around 1360 MeV/c^ 
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16l | and which is usually buried by the strong A peak in vrA^ and experiments; 



the other one around 2065 MeV/c^ may be one or more of the long sought "missing" A^* 
resonance(s). 

To estimate the mass and width of these resonance peaks, we use simple relativistic Breit- 
Wigner functions plus a smooth function to fit the data points of Fig. El For the final fit we 
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FIG. 6: Data divided by Monte Carlo phase space versus pir invariant mass for J/ip ^ pir^n (solid 
circle) and J/tlj — > pvr^n (open square), compared with our fit (solid curve). The contributions of 
each resonance peak are shown by the dot-dashed lines in the same figure. The dashed line is the 
contribution of background terms including the nucleon pole term. 

use 

^ I ^ I 

\A\' = |Co + C',s^,\ + g (,^^ _ Mfy + MfVj 

where s.^-^ = M^^ is the invariant mass squared of the pvr system. Background is described 
with a constant and s-i^p linear term plus an additional Breit-Wigner at low energy to simulate 
the contribution from the nucleon pole term and the A and other backgrounds. All param- 
eters are free. The fitted mass and width for the four A^* peaks are listed in Table HI The 
systematic error in the table is obtained from the variation of fit results when assuming other 
background functions, such as |Co|, \CiMt,p\, \C'qS.„p\, \Cq + CiM^rpI, \Cq + CiM^^p + CqS^^-pI, 
and adding one more A^* resonance around 1.85 GeV/c^ or above 2.2 GeV/c^ in the fitting. 
The fit to the data is also shown in Fig. IHlwith = 133 for 102 data points. All four peaks 
are highly significant. The statistical significance is 11a for the A^*(1440) peak (the least 
significant one) and 13cr for the new A^*(2065) peak. 
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TABLE I: The fitted masses and widths for the four A^* peaks shown in Fig. |H1 



Mass(MeV/c2) 


Width(MeV/c2) 


1358 ± 6 ± 16 
1495 ± 2 ± 3 
1674 ± 3 ± 4 
2068 ± 3tll 


179 ± 26 ± 50 
87± 7± 10 
100 ± 9 ± 15 
165 ± 14 ± 40 



Because we use a constant width for the Breit-Wigner (BW) formulae, the BW mass and 
width are very close to their corresponding pole positions. For the two well-known peaks 
at 1500 MeV/c^ and 1670 MeV/c^, the former contains two well-established N* resonances, 
the iV*(1520) and A^*(1535), while the latter contains more than two A^* resonances. Here 
we only use one BW function to fit each of them. 

For the new A^*(2065) peak, orbital angular momentum L = is preferred due to the 
suppression of the centrifugal barrier factor for L >1. If we fit the A^*(2065) peak in Fig. IHl 
with L = 1 centrifugal barrier factor instead of Eq.(l), then the increases from 133 to 
163 for 102 data points. The much worse fit with L = 1 compared with L = means that 
there is substantial L = component for the new A^*(2065) peak. 

For L = 0, the spin-parity of A^*(2065) is limited to be 1/2"*" and 3/2+. This may be 
the reason that the A^*(2065) shows up as a peak in J/ip decays while no peak shows up in 
the TiN invariant mass spectra in vrA^ and 7A^ production processes which allow all 1/2^, 
3/2^, 5/2^, and 7/2^ A^* resonances and their isospin 3/2 A* partners around this energy to 
interfere with each other. In order to determine its spin-parity, a partial wave analysis using 
an effective Lagrangian approach is tried for the pTT~n data by including a new A^*(2065) 
with spin-parity either 1/2+ or 3/2"*" in addition to all well-established A^* resonances below 
2 GeV/c^ with masses and widths fixed to their PDG values Q]. Comparing with the fit 
without including any new A^*(2065) resonance, including a A^*(2065) with spin-parity of 
either 1/2+ or 3/2+ improves log likelihood value by more than 400. The new A^*(2065) 
peak cannot be reproduced by reflections of well-established A^* resonances. However, the 
spin-parity of the new resonance(s) cannot be well determined. The difference of the fits 
with different spin-parity quantum numbers is small and depends on many fitting details, 
such as how to treat the background contribution and how large an isospin breaking effect is 
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allowed. Including both 1/2"*" and 3/2"'" improves log likelihood value further by more than 
60. So it is quite possible that both are needed. There are quark model predictions for the 
existence of A^* resonances with spin-parity 1/2"'" and 3/2"'" in the energy ran ge f rom 2.0 to 
2.1 GeV/c2 0,0,0. The same channel has also been studied in ip' decays |19|. However 
due to lower statistics and the allowance of more partial waves, only a much broader peak 
was observed above 2 GeV/c^. 

In summary, the J/tp ^ NNn decay at BEPC provides an excellent place for study- 
ing pure isospin 1/2 excited A^* resonances. Using 58 million J/ip decays, more than 100 
thousand J/ip ^ pix^n + c.c. events are obtained. The corresponding branching ratios are 
determined to be (2.36 ± 0.02 ± 0.21) x lO^^ and (2.47 ± 0.02 ± 0.24) x 10"^ for p-n-fi and 
p7T~^n, respectively. In the pn invariant mass, besides the two well-known peaks around 1500 
MeV/c^ and 1670 MeV/c^, there are two new clear peaks. The one around 1360 MeV/c^ 
with statistical significance of llo", identified as due to the A^*(1440) resonance (usually 
obscured by the strong A peak in ttN and experiments), has a Breit-Wigner mass and 
width of 1358 ± 6 ± 16 MeV/c^ and 179 ± 26 ± 50 MeV/c^. The other N* peak around 2065 
MeV/c^ with statistical significance of 13cr is believed to be due to long sought "missing" 
A^* resonance (s) predicted by many theoretical quark models in this energy range. A 
simple Breit-Wigner fit gives its mass as 2068 ± 3!l^4o MeV/c^ and width as 165 ± 14 ± 40 
MeV/c^. From a partial wave analysis, we conclude that the new A^*(2065) peak cannot be 
reproduced by reflections of PDG established A^* resonances and includes at least one new 
A^* resonance. 
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